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Abstract 

The paper addresses the problem of a Mode III interfacial crack advancing quasi-statically in a 
heterogeneous composite material, that is a two-phase material containing elastic inclusions, both 
soft and stiff, and defects, such as microcracks, rigid line inclusions and voids. It is assumed that 
the bonding between dissimilar elastic materials is weak so that the interface is a preferential path 
for the crack. The perturbation analysis is made possible by means of the fundamental solutions 
(symmetric and skew-symmetric weight functions) derived in Piccolroaz et al. (2009). We derive the 
dipole matrices of the defects in question and use the corresponding dipole fields to evaluate effective 
tractions along the crack faces and interface to describe the interaction between the main interfacial 
crack and the defects. For a stable propagation of the crack, the perturbation of the stress intensity 
factor induced by the defects is then balanced by the elongation of the crack along the interface, 
thus giving an explicit asymptotic formula for the calculation of the crack advance. The method 
is general and applicable to interfacial cracks with general distributed loading on the crack faces, 
taking into account possible asymmetry in the boundary conditions. 

The analytical results are used to analyse the shielding and amplification effects of various types of 
defects in different configurations. Numerical computations based on the explicit analytical formulae 
allows for the analysis of crack propagation and arrest. 
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1 Introduction 



Analytic solutions for a crack propagating in a homogeneous elastic solid containing a finite number of 
small defects (elastic and rigid inclusions, microcracks, voids) have been derived in Bigoni et al. (1998), 
Valentini et al. (1999), Movchan et al. (2002) on the basis of the dipole matrix and weight function 
approach. The dipole field describes the leading-order perturbation produced by a small defect placed 
in a smooth stress field and gives rise to "effective" tractions applied along the crack faces. Furthermore, 
the weight functions allow for the derivation of the corresponding perturbation of the stress intensity 
factor as weighted integral of the "effective" loading. The method is general and applicable to both 2D 
and 3D cases and to defects of different type and shape, provided that the corresponding dipole matrix 
is appropriately constructed and the weight functions for the corresponding unperturbed cracked body 
are available. 

Problems of a macrocrack interacting with microcracks have been analysed by Romalis and Tamuzh 
(1984) under the infiuence of mechanical loading and by Tamusz et al. (1993) under the influence of heat 
fiux, using the analytic functions and singular integrals approach (Muskhelishvili, 2008). The possible 
closure of crack surfaces and the consequent appearance of a contact zone have been considered in 
Tamuzs et al. (1994) and Tamuzs et al. (1996). Asymptotic models of a semi-infinite crack interacting 
with microcracks have been developed by Gong and Horii (1989) and Meguid et al. (1991) using complex 
potentials and the superposition principle. A review of publications on macro-microcrack interaction 
problem is given in Tamuzs and Petrova (2002). 

Recently, Piccolroaz et al. (2009) derived the symmetric and skew-symmetric weight functions for 
a semi-infinite two-dimensional interfacial crack, thus disclosing the possibility of applying the dipole 
matrix approach to the propagation of cracks along the interface in heterogeneous materials with small 
defects. The weight functions constructed in Piccolroaz et al. (2009) are of the generalised type and 
thus applicable to any type of boundary conditions along the crack faces. The singular perturbation 
associated with a small crack advance was also obtained there. It is worth noting that the availability of 
the skew-symmetric function for the problem under consideration is essential since the "effective" loading 
produced by the defects on the crack faces is not symmetrical, in general. In the present paper, we analyse 
the scalar case of antiplane shear loading. The full vector problem will be addressed elsewhere. 

The paper is organised as follows. The formulation of the problem is outlined in Section [21 which 
includes also preliminary results on the unperturbed problem. Section |3] is devoted to the perturbation 
analysis, in particular the derivation of the dipole fields for several types of defects and the analysis 
of singular perturbation associated with the crack advance. Section |4] provides a number of numerical 
results based on the explicit analytical formulae. In particular, shielding and amplification effects of 
different defect configurations on the crack-tip field are presented. The possibility to design a neutral 
configuration for any given force system distributed along the crack faces is discussed. Finally, the crack 
propagation and arrest produced by the defects under consideration is analysed. In the appendix, we 
derive the dipole matrix for an elliptic elastic inclusion placed in a homogeneous antiplane field. 

2 Problem formulation and preliminary results 

We consider a two-dimensional composite structure consisting of a bimaterial matrix (two dissimilar 
elastic half-planes il±) containing a dilute distribution of inclusions, microcracks and rigid line inclusions, 
see Fig.[T] The two materials constituting the matrix are assumed to be linear elastic and isotropic, with 
shear moduli denoted by and respectively. All interfaces between different phases are assumed 
to be perfect, that is, the displacements and tractions remain continuous across the interface. 

We introduce the following notations. Let C r2+ be a small elastic inclusion of diameter 2eli 
centred at the point Yi ~ (ai, hi). The shear modulus of the inclusion is denoted by /^i, and it can be 
greater or smaller than the shear modulus fi-^- of the surrounding material, so that both stiff and soft 
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inclusions are considered. The notation 7! C fi- is used for a microcrack of the length 2el2, centred 
at the point Y2 = (02, /12) and making an angle a2 > with the positive direction of the a;i-axis. By 
7I C f2+ we denote a movable rigid line inclusion of the length 2elz, centred at the point = (03, /is) 
and making an angle as > with the positive direction of the xi-axis. Although these notations refer 
specifically to Fig. [1] the formulation can be easily extended to problems with different number and type 
of defects, as those considered in Sec. 21 

We assume that a semi-infinite interface crack advances quasi-statically along the interface 
connecting the half-planes, and we denote the uniform advance of the crack by e^cf). Here and in the 
sequel, e > is a small dimensionless parameter. The reason for the order in the crack advance will 
be clear in Sec. 121 where we perform the asymptotic analysis. 
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Figure 1: Geometry of the problem: interface crack in a bimaterial plane with defects: denotes a small elastic inclusion, 
7| a microcrack, and 7| a rigid line inclusion; Yi, Y2, ^3 are the "centres" of the defects. 



We assume that the composite is dilute, that is, the small defects are distant from each other so 
that the interaction between them can be neglected. Consequently, we can model the three cases of an 
elastic inclusion, of a microcrack and of a rigid line inclusion separately. It is possible then to extend the 
results to a finite number of defects by superposition using the linearity of the problem, provided that 
the distance between defects remains finite. 

An external loading p± is applied to the crack faces Fj_ and it is assumed to be self-balanced such 
that the principal force vector is zero, that is, 

p+dx — / p^dx = 0. (1) 
Jrt 

We assume that the loading p± on the crack faces vanishes in a neighbourhood of the crack tip. 
The problem is then formulated in terms of the Laplace equation 

AM±(xi,a;2) = 0, Aui{xi,X2) = 0, (2) 

where u = {u^,u_,Ui} denotes the displacement component along xs-axis in the respective domain 
fl+ \ {g^ U 7|}, n_ \ 7f and g,. 

We prescribe the following boundary conditions on the crack faces 

fJ'±^ — =P± onF±, (3) 
0x2 
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and the interface conditions 



u+ = u-, n+- 



on 



dx2 dx2 

The transmission conditions for the elastic inclusion are formulated similarly to that is, 

du+ dui 
u+=Ui, ^1+^^ = ^ii— on age- 

We assume that the microcrack faces 7^ are traction free, so that 



9u_ 
dn 



= on 72 



(4) 



(5) 



(6) 



Finally, on the boundary of the movable rigid line inclusion 7I the Dirichlet boundary condition is 
prescribed, that is, 

u+ = on 7I, (7) 
where is an unknown constant which will be defined later from the balance condition 



du+ 
dn 



ds ^ 0. 



(8) 



The unperturbed problem (e = 0) corresponds to a semi-infinite interfacial crack in a bimaterial 
plane. The solution to this problem (see, for example, Piccolroaz et al. 2010) can be expressed in the 
form 



uf{r.e) = — / u^l\sSr-'ds, ±6 G (0,7r), 



where 



sin(s6') 
/i±s cos(7rs) 



cos(.s6 



{fi+ — fi-) sin(.s6 



(/i+ + /i-)s sin(7rs) 2/i±(/i+ + /i_)s cos(7rs) 



(9) 



(10) 



and < Lo < 0.5, with the symmetric and skew-symmetric components of the loading denoted by 
W ^ P+ ^ P- and (p) = (p+ + p-)/2, respectively. 

The two-terms asymptotic expansions for the traction a^*^'^ ahead of the crack tip and for the crack 
opening lu^^^ read as follows 
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/ {(p)(-r) + |bK-r)}r-i/2dr, 
I r {{p)i-r) + lMi-r)}r-^/'dr, 



(11) 
(12) 

(13) 
(14) 



with 77 = (/_*_ — i^i^)/{fi^ + /i_) being the contrast parameter. 
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For the perturbation analysis we will also need the so-called symmetric and skew- symmetric weight 
functions, [[f/] and ([/), respectively, and the corresponding traction (E) in a bimaterial plane with a 
semi-infinite crack (see Piccolroaz et al., 2009), 



f - i 



-1/2 



, xi > 0, 
xi < 0, 



xi > 0, 



(1 - i)fJ-+^J- 



(-a-i)-3/2, XI < 0. 



(15) 



(16) 



2V27r(/i+ +^_) 

For an arbitrary solution u of the governing equation ^ and the corresponding traction a along the 
xi-axis, the reciprocity identity has the form (Piccolroaz et al., 2009) 

{lUiix[ - x,){a)ix,) + {U)ix[ - Xi)H(a^i) - (S)(a;; - x,)Mix,)}dx, = 0, (17) 

where the following notations have been used: 

[/! = /+-/_, (/) = !(/+ + /_). 

Let us introduce the notations 

/(+)(xi) = /(a;i)F(xi), f^-\xi) = f{x^)H{-x^), 
where H denotes the Heaviside function, so that 

/(a;i) = /(+)(xi) + /(-)(xi). 
The reciprocity identity (jl7p can be written as 

{[C/l(x; - x^){a)^+Hx,) - (S)(.t; - x,)luf -\x^)} dxr = 

{ iC/Kx'i - xi)(a)(-)(xi) + (C/)(x'i - x,)laf-\x,)\dx^ 

oo 

{U){x[ - xi)l<jf+Hxi)dxi + / (S])(x'i - xi)luf+\xi)dxi. 



(18) 



(19) 



Note that (cr)^^-', Icr]^^-* are the average and the jump of the prescribed loading on the crack faces, 
whereas [cl lu]^^-* are the prescribed discontinuities of traction and displacement across the interface. 
The identity ([T9|) will be used extensively later in the perturbation analysis. 



Additionally, we shall need the components of the displacement gradient Va;w^'' computed at an 
arbitrary point Y = {dcosip, dsimp). These are obtained from and given by 
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dx2 



Y 



1 



dxi 



2 cos Lp — x\l d — dj x\ 
, 2(p)(xi)+7/ 




(21) 



3 Perturbation analysis 

We shall construct an asymptotic solution of the problem using the method of Movchan and Movchan 
(1995), that is the asymptotics of the solution will be taken in the form 

3 3 

u±(a;,£) =ui"^(a;) + £^W^j(4j)+£^^ti2'^(a;)+£^w(a;,0)+o(£2), e ^ o. (22) 

In (im, the leading term (x) corresponds to the unperturbed solution, and it is described in the 
previous section. The term £X]^=i ^j(^j) corresponds to the boundary layers concentrated near the 
defects and needed to satisfy the respective conditions ([5]), © and ([7|); the variables will be defined 

in the next section. The term u'± {^') is introduced to fulfil the original boundary conditions 

on the crack faces and the interface conditions (j?]) disturbed by the boundary layers; this term, in turn, 
will produce perturbations of the crack tip fields and correspondingly of the stress intensity factor. In 
Sec. 13.11 we will analyse the effect of each defect separately. Finally, the term e^vix, (p) corresponds to a 
singular perturbation associated with a possible crack advance, £^0, and will be addressed in Sec. 13.21 



3.1 Singular perturbation and dipole fields generated by small defects 
3.1.1 Small elastic inclusion 

We shall start with the elastic inclusion, situated in the upper half-plane. The leading term 

clearly does not satisfy the transmission conditions ([5|) on the boundary dg^- Thus, we shall correct the 

solution by constructing the boundary layer Vl^i(^i), where the new scaled variable is defined by 



(23) 



with Yi being the "centre" of the inclusion t/g (see Fig. [Ij. 

For 1^1 (li) = {W^i™, e g; W^"\ |i S \ g} we consider the following problem 

AI^r(^i)-0, liGg, AWr'^O, ^,eR'\g, (24) 

where 

g = e-'g, = e . ^ ^ 
The function Wi remains continuous across the interface dg, that is, 

Wl" = on dg, 

and satisfies on dg the following transmission condition 

t^,-^Wr{i,)-fi+-^WrHii) = {f^+-f^.)n-Vu^'HYi) + Ois), £^0, (25) 
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where n = is an outward unit normal on dg. The formulation is completed by setting the following 
condition at infinity 

M^°"*^0 as llil^oo. (26) 

The problem above has been solved by various techniques and the solution can be found, for example, 
in Movchan and Movchan (1995). 

Since we assume that the inclusion is at a finite distance from the interface between the half-planes, 
we shall only need the leading term of the asymptotics of the solution at infinity. This term reads as 
follows 



1 



(0) 



1^1 



Ml 



^1 



l^il 



+ 0(|li|-2) as li^oo. 



(27) 



where A4i is a 2 x 2 matrix which depends on the characteristic size li of the domain g and the ratio 
12+/ Hi] it is called the dipole matrix. For example, in the case of an elliptic inclusion with the semi-axes 
la and If, making an angle ai with the positive direction of the xi-axis and a;2-axis, respectively, the 
matrix A^i takes the form 



Mi = --lalb{l + e){fi.-l} 



1 + cos2q;i 1— cos2q:i (1 — e)(/ii, — 1) sin2ai 
e + /-i* 1 + e^* {e + fii,){l + e^J.^,) 

(1 — e)(/i* — 1) sin 2q;i 1 — cos2ai 1 + cos2q!i 



1 



(28) 



where e = h/la and /x+ = fi+/fii. We note that for a soft inclusion, ^+ > /ii, the dipole matrix is 
negative definite, whereas for a stiff inclusion, < /li, the dipole matrix is positive definite. 

The term eWi{^i) in a neighbourhood of the .Ti-axis written in the x coordinates takes the form 



(29) 



where 



Ml- 



X 



Yi 



Yi\ 



(30) 



As a result, one can compute the average e^(<Ti) and the jump e'^|cri] of the "effective" traction across 
the line X2 = induced by the elastic inclusion g^. Since du^^'^ /dx2 = —dwi/dx2 must hold on the line 
2^2 = (to satisfy the original boundary conditions ^ and interface conditions (U)), this gives 



(-l>=-2(M++M-)5^: 



dwi 



(31) 



where 



dwi 
dx2 



1 

2^ 



Yr 



Ml 



62 



Y^ 



In the limit fii — > cxj, we obtain the dipole matrix for a rigid movable inclusion. In the case of an 
elliptic rigid inclusion, we have 



M 



-Ubil/e + 1) 



1 + cos 2a + e(l — cos 2a) (1 — e) sin 2a 

(1 — e) sin 2a 1 — cos 2q! + e(l + cos 2a) 



(33) 



8 



3.1.2 Microcrack or small void 



Now we apply the same procedure to a microcrack. Since the leading term u^^") does not satisfy traction 
free boundary conditions © on the microcrack faces 7^, we construct a boundary layer 1^2(^2)1 where 
the new scaled variable ^2 is defined similarly to with Yi replaced by Y2, the middle point of the 
microcrack (see Fig. [T]). 

The function W2 satisfies the Laplace equation in a homogeneous plane containing a crack of finite 
length 2I2 and the following traction conditions on the crack faces 



-n-Vu^°\Y2) + 0{e 



0, 



(34) 



where n 
: 



is an outward unit normal on the crack faces of the finite crack 72 = £ 72 = {^2 ^ 

Y2 €72}- Again, we are looking for the boundary layer solution which decays at infinity. 
Such a problem has been solved by various techniques and the solution can be found, for example, in 
Muskhelishvili (2008). The asymptotic formulae ((27|) . ([29| and ((30|) remain the same with the dipole 
matrix Aii replaced by 



M2 



sm a2 
— sin a2 cos a2 cos^ a2 



— sm Q!2 cosq;2 



nl2 



1 — cos 2a2 — sin 2a2 
— sin 2^2 1 + cos 2q;2 



(35) 



giving the average (cr2) and the jump [[(T2I of traction across the interface between the two half-planes 
in the form of pip with wi replaced by W2- 

For the case of a small void, one can construct the corresponding boundary layer following Movchan 
et al. (2002) . It has been shown that for an arbitrary finite void, the dipole matrix coincide with the one 
corresponding to an "equivalent" elliptic void. Denoting by la and 4 the semi-axes of the ellipse, and by 
a2 the angle formed by ^a-semi-axis and the xi-axis, the dipole matrix can be written as follows 



M 



-ila+lbf 



1 - e 



cos 2a — sin 2a 
sin 2a cos2q! 



where e = lb/la '^s the eccentricity of the ellipse. One can see that in the limit case as e 
(j35p is recovered from 



(36) 



0, the formula 



3.1.3 Rigid line inclusion 

In the case of a small rigid line inclusion, the boundary layer solution 1^3(^3) satisfies the Laplace 
equation in the entire plane of the scaled coordinates system ^3 G introduced as in ([2^. where 
the point Yi is replaced by 1^3 and stands for the middle point of the rigid inclusion (see Fig. [T|). 

On the boundary of the rigid inclusion 73 = £"^73 = {^3 G : £^3 + 1^3 € 73} the function W3 
satisfies the following condition (compare with ([7])): 



£^3(13) = - u^^Hy^) - e ^xu^"^ 
Substituting 1^ into © gives 



l3 + 0(£2), £^0. 



= ufiY,)+eVxU^+\, •C3 + £W^3(l3)+0(£2), £^0. 



To satisfy this condition we need to set u, = u^_^\Y3) and 



^3(^3) = - Vxu 



(0) 



C3> ^3 € 73- 



(37) 



(38) 
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Finally, we are looking for a solution which decays at infinity similarly to (P^ . 

Such a problem can been solved by various techniques. We shall present the construction of the 
dipole field using complex analysis in [B] following Muskhelishvili (2008). The corresponding dipole 
matrix takes the form 

1 + cos 2a3 sin 2a3 



sin 2a3 



1 — cos 2a,-! 



(39) 



One can see that in the limit case as e — >■ 0, the formula (p9| is recovered from ([33)1 . 

In the original coordinate system x, far away from the rigid inclusion, the asymptotics of the boundary 
layer solution W3 takes the same form as (PH)) and (PU)) . where wi and Aii should be replaced by W3 
and M-Z; respectively. 

The average and jump of traction across the interface crack can be found by formulae analogous to 



3.1.4 Line defect with imperfect bonding 

The dipole field and corresponding dipole matrix for a line defect with imperfect bonding can be derived, 
as a limiting case, from the solution for a perfectly bonded elastic elliptic inclusion, see Sec. 13.1.11 We 
distinguish between stiff and soft line defect. 

The transmission conditions for a soft line defect (Antipov et al.. 2001; Mishuris. 2001; Hashin. 2001) 
are given by 

H(s)=0, lul{s) = Ka[s), (40) 

where s is an abscissa along the line inclusion, \u\ and a are the displacement jump across the line inclu- 
sion and the respective (continuous) traction, and k, is the compliance of the bonding. Correspondingly, 
the dipole matrix is derived from (j28[) by choosing fi^lf, — k and then taking the limit If, ^ 0. We obtain 



M 



21 



1 



- cos 2a 
sin 2a 



1 



sin 2a 
I- cos 2a 



(41) 



One can check that the solution for a microcrack psp is recovered by taking the limit k cx) in (|¥T|l . 

For a stiff line defect, the transmission conditions (Benveniste and Miloh, 2001; Mishuris, 2003) are 
written as 



0. 



In this case, the dipole matrix is derived from 
lb — >■ 0. We obtain 



by choosing /i* 



I' 



2 1 + kI 



1 + cos 2a 
sin 2a 



sin 2a 
1 — cos 2a 



(42) 

nit and then taking the limit 

(43) 



One can check that the solution for a rigid line inclusion ppp is recovered by taking the limit k — > in 

(1131). 



3.2 Singular perturbation for crack advance 

In this section, we consider a singular perturbation of the physical fields generated by an advance of 
the crack tip by a small quantity e'^ip. We denote unperturbed quantities by subscript and perturbed 
quantities by subscript We also write the perturbed fields with reference to a coordinate system 
centred at the new crack tip position. 
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Writing the Betti identity (fTO)) for the unperturbed ((Tq and |mo]) and the perturbed (cr^ and Jw*]]) 
quantites and then subtracting one from the other, we obtain, after the Fourier transform, (see for details 
Piccolroaz et al., 2009) 

[C71 + K+ - e'^^' V) - S"(Nr - e^^'^'^m-) - 0. (44) 

In the bar denotes Fourier transform with respect to the variable xi, and /3 the corresponding 
variable in the transformed space. The superscripts + and — characterise functions analytic in the 
upper and lower half planes respectively. 

Following Willis and Movchan (1995), we expand the exponential term as exp(i/3e^0) = 1 + ie'^(f)(5 + 
0(e*) and also substitute into (|44p the two-terms asymptotic of traction and crack opening [[uo]^, 
as /3± — > oo: 

^ _ (l + ^)(A^+ + /^-)< ^-3/2 ^ (l-^)(AH+M-)<^ ^-5/2 ^ ^ r/2 
2^+^_ 4/i+/^_ 

Note that the perturbed fields, ct+ and [[u*!^, have the same asymptotic expansion as in (|45)) and (|46)) . 
subject to replacing _ftr|J^j with and with Ajjj. Then, collecting like powers of /3±, we get 



2 



(^4°^ - /^m) + ^-J^^e'ct.A\,^ {(3^' pz') + 0{p-') = 0. (47) 



From P7)l we can conclude that, in the limit as e — > 0, the first-order perturbation of the stress intensity 
factor due to a uniform advance of the crack tip is given by 



i^f„-<':=e^A< = i/<\ (48) 



where ^jjj is the coefficient in the second-order term asymptotics of the unperturbed fields and it is 
given in terms of the loading applied to the crack faces as follows 







<i ^ + iM(^i)) i-^ir'^'dx,. (49) 

Note that, although a competing double expansion, as e — !• and (3± — >■ oo, has been used in this 
procedure, the correctness of the final result has been proved in Piccolroaz et al. (2009) by means of the 
Wiener-Hopf technique. 

3.3 Analysis of a stable quasi-static propagation of an interfacial crack 

The stress intensity factor is expanded as follows 

Km = <i + I^AXfjj + J2 ^^m j + o{e'), e 0, (50) 

where k'^^ is the stress intensity factor for the unperturbed crack, Aifj'^j is the singular perturbation 
due to the uniform advance of the crack tip, ^-^m perturbation produced by the defects. 
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Assuming that the crack is advancing quasi-statically along the interface, the energy release rate 
remains constant and equal to the critical value, G = Gc, so that 



AG = G - G(°) = 0. 

The energy release rate can be represented in terms of the stress intensity factor by 

G 

so that ([?T|) becomes 



-(- 

4 V/i+ 



1 



R'2 



A' 



(0) 



^III - -'Mil ■ 

Upon replacing in the above equation the expansion (j50p . we obtain 



(51) 

(52) 
(53) 



2£2<|A< + ^AA-(^^|+o(e^)=0, 



and thus 



AK 



III 



where Ai^j^Jj is given by (|35]), whereas AA'|jj is given by 



U) _ 
III ~ 



0, 



(54) 



(55) 



(56) 



with (cr'-''^) and fcr^-'^ being the average and the jump of the traction induced by the defects across 
the crack faces, given by ([5T |) - ([5^ . The integral in ([55]) can be evaluated explicitly and we obtain the 
formula 



AK 



III 



2 pt+Ai- 



(57) 



where M.j is the dipole matrix for the corresponding defect, Vaju'^^jy^ is the gradient of the displace- 
ment in the unperturbed configuration evaluated at the centre of the defect Y j — {dj cos (pj,dj simpj) 
given by (PT|) . and 



2df 



- sm — ^, cos — r- 
2 ' 2 



(58) 



As a result, we obtain the formula 



A 



(0) 

III 3 = 1 



III ' 



(59) 



which, together with (1^5)1 and ([57)1 , allows for the analysis of the quasi-static propagation of the interfacial 
crack interacting with small micro-cracks and inclusions. 



12 



4 Numerical results 



In this section, we show numerical resuUs based on the exphcit asymptotic formulae for the perturbation 
of stress intensity factor ([57)) and for the crack advance (|59p . In Sec. 14.11 we analyse the shielding and 
amplification effects on the crack-tip field of two defect arrangements (see Fig. [5]) . From this analysis 
it is possible to predict whether the crack will propagate (amplification) or not (shielding) from the 
specified position. In the case of propagation, the formula (|59p allows us to incrementally calculate the 
crack advance, so that we can predict the extent of crack propagation and, in particular, whether the 
crack will arrest at a certain point or continue propagating. This is discussed in Sec. 14.21 



4.1 Shielding and amplification eff^ects. 

A neutral defect configuration is a configuration for which the perturbation of the stress intensity factor 
is zero, AKm = 0. For non-neutral configurations, we may have shielding effect when the defects produce 
a decrease of the stress intensity factor, AA'ni/A'|^j < 0, or amplification effect in the opposite case, 
AKui/kI^} > 0. 

Before the discussion of computations, we derive simplified asymptotic formulae, which are valid in 
the case where the loading is applied far away from the crack tip. Let us denote the support of the 
applied loading by supp± C (— cxi,— a), a > 0. We have the following reduced asymptotic formula as 
a — >■ +00, valid for any elliptic inclusion with semi-axes la,lb and making an angle a with the positive 
a; 1 -direction 



A/\ni llJb,^ , w /^T / 1 • /3v3 \ ■ 



"^iii 



1 f3ip 



1 + e/i* V 2 



(60) 



From this formula we can derive the limiting cases as a — > +oo for different type of defect, namely: 

- for a microcrack: 

AA'iii 1 l^ A^T /^3<^ \ \ /flix 

TTTT TT-TT COS — Q; COS — — Q! : (61) 

- for a rigid line inclusion: 
AAm., IP .„^3^ \. (62) 



2dV++/^- V2 ; V2 

- for an elliptic void: 

^ ~ - "?7^ - ") - »' - ") » (f - ") } ^ 

- for an elliptic rigid inclusion: 

AA'iii ILk,^, , ix / ■ /^^(^ \ ■ fV' \ , f^f \ 

TTTT- ~ — 1 / e + 1 ' < sm a sm a] + e cos a cos a 

K^o) 2 rf2 W V2 J V2 J \2 J \2 

(64) 

- for a line defect with soft bonding and transmission conditions given by (HO]): 



AA'iii 1 r u=r k (Zw 

■ cos — a I cos 



k'^) 2d^ IJ,+ + li-l + K V 2 
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- for a stiff line defect with transmission conditions given by 



2 i^L+ + fi- 1 + kI V 2 / V 2 

Note that the formula in the case of a homogeneous body was previously obtained by Gong (1995). 

From these simplified asymptotic formulae, it is possible to design neutral configurations, for the case 
where the loading is applied at large distance from the crack tip. To illustrate this, we show in Fig. [2] 
two examples of defects arrangements. The case of Fig. [5^ consists of two defects, a microcrack, denoted 
by subscript 1, and a rigid line inclusion, denoted by subscript 2, located in the same half-plane and 
such that 

h/di^l2/d2, 01-02=0, ai-a2 = Tr/2. (67) 

The case of Fig. ^jp consists of the same two defects as above, a microcrack and a rigid line inclusion, 
but they are located in different half-planes and such that 

H-ll/dl^ H+ll/dl, 01+02 = 0, ai+a2=7r/2. (68) 




Figure 2: Two neutral arrangements for remotely applied loading, a — > +oo. Arrangement (a): a microcrack of length 
2£l\ , orientation oi and placed at the point (di cos ipi , d\ sini/Ji) and a rigid line inclusion of length 2ei2 = 4£ii , orientation 
02 = ot\ — 7r/2 and placed at the point {2d\ cos ip\ , 2d\ sin i^i). Arrangement (b): a microcrack of length 2£l\ , orientation 
oil and placed at the point (di cos y'l , di sini^i) and a rigid line inclusion of length 2el2 = 2eli ^ /i- , orientation 
Oil = t/2 — oi and placed at the point (di cosi/3i, — di sin^Pi). 



When the loading is applied at a finite distance from the crack tip, arrangements (a) and (b) shown 
in Fig. [5] are, in general, not neutral any longer, and we may have shielding and amplification effects. 
To illustrate this, we consider a three-point load as shown in Fig. (5] This load is self-balanced for any 
b < a, and is symmetrical only for 6 = 0. 

We consider first arrangement (a) for the case of a weak interface between two identical materials 
(/i_i- = = yu) and symmetrical load, 6 = 0, see Fig. [31 In the computation we use the values: sli = 0.1, 
di = 1, e/2 = 0.2, d2 = 2. On the diagrams, the horizontal axis stands for the angle ipi (— tt < ipi < tt) 
defining the angular position of the centre of the microcrack 1 with respect to the crack tip. On the 
vertical axis we report the value of the angle ai (0 < ai < tt) defining the crack orientation with respect 
to the xi-direction. The position of the rigid line inclusion 2 is determined by ()67|) . It is worth noting 
that, in a real experiment, the perturbation of the stress intensity factor AKm can be measured with the 
accuracy 6 = A/ini/A'|n with respect to the unperturbed value Kjjj, so that the region corresponding 
to shielding is defined by AKm < ^SK^j^ and shadowed by light grey in the diagrams. The region 
corresponding to amplification is defined by AiiTiii > 5K^^^ and shadowed by medium grey. The dark 
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grey region in the diagrams corresponds to situations for which |A_ftrni| < '^^m ; that is the perturbation 
is neghgibly small and the corresponding configuration could be regarded as neutral. 

Fig. [3K corresponds to a two-point (symmetrical) load situated close to the crack tip, a — 3. The 
diagram is symmetrical with respect to the angle ipi, as the load is symmetrical and the material is 
homogeneous. As the distance a of the load from the crack tip increases, the diagram changes and the 
dark grey region, corresponding to neutral configurations, enlarges, as shown in Fig. [SJj where we set 
a = 100. 




(a) (b) 



Figure 3: Shielding (light grey), amplification (medium grey) and neutral (dark grey) regions created by a microcrack 
and a rigid line inclusion arranged as shown in Fig. [2ji. The diagrams are for a weak interface between two identical 
materials, rj = 0, and symmetrical loading, b = 0. The position of the microcrack is characterised by the angles ipi and 
oi, whereas the position of the rigid line inclusion is determined by conditions 1671 . The dark grey region corresponds to 
neutral configurations, that is configurations for which the perturbation is negligibly small (within the accuracy S = 10~^). 
Figure (a): load is applied close to the crack tip, a = 3. Figure (b): load is applied at a large distance from the crack tip, 
a = 100. 

In Fig. 01 we analyse the effect of the asymmetry of loading for the same configuration as in Fig. |31 
The diagrams correspond to a three-point loading, characterised by the parameters a and b, as shown 
in Fig. [2l The definition and range of the angles ai and ipi are the same as in Fig. [3] (— tt < ipi < it, 
< ai < tt). We consider four values of the parameter a, namely a = 3, 6, 10, 100, from left to right. This 
allows us to trace different situations corresponding to the increasing distance of the loading from the 
crack tip (keeping b fixed) . The asymmetry of loading is measured by the parameter 6, which increases 
from the upper part (symmetrical loading 5 = 0) to the lower part (5 = 0.99a) of the figure. 

From the results presented in Fig.|4l the following conclusions can be drawn. First, the neutral region 
(shadowed dark grey) enlarges when moving from the diagrams on the left to those on the right. This 
finding, already evident in Fig. [3] where only the values a = 3 and a ~ 100 were considered with 6 = 0, 
confirms the validity of the simplified formulae (j6l]) and (|62p . established in the literature (Gong, 1995). 
Secondly, it is clear that the asymmetry of loading plays an important role. Indeed, even for a = 100 
and b = 0.33a, or 6 = 0.66a, when the nearest point force is still well separated from the crack tip, the 
shielding-amplification diagrams are essentially different from those of the symmetrical case (6 = 0). 
Thirdly, since the distance of the loading from the crack tip is measured by a — 5, the neutral region 
shrinks by increasing the asymmetry b at fixed a (thus moving from the diagrams on the upper part to 
those on the lower part). 

In Fig. [S] and Fig. [SI we show the infiuence of the inhomogeneity on the shielding and amplification 
effects. For this purpose, we consider the aforementioned three-point loading for the case of a bimaterial 
plane with contrast parameter 77 = 0.67 (diagrams of Fig. [5]) and 77 = —0.67 (diagrams of Fig. [6]). All 
other parameters (a, 6, ai and ipi) are defined in the same manner as above. 

The infiuence of the inhomogeneity is clearly observable only for a highly pronounced non-symmetrical 
loading which is applied close to the crack tip. Furthermore, the shielding and amplification effects are 
more pronounced when the defects are located in the softer material (upper half-plane, < ipi < tt, in 
Fig.m lower half-plane, — tt < t^i < 0, in Fig.[6|). 

Finally, we comment on the arrangement (b), see Fig. [2)d. This configuration turns out to be neutral 
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a = 3 fl = 6 a =10 «=100 




V>\ t-l tfl Wl 

Figure 4: Shielding (light grey), amplification (medium grey) and neutral (dark grey) regions created by a microcrack 
and a rigid line inclusion arranged as shown in Fig. [2^. The diagrams are for a weak interface between two identical 
materials, ?J = 0, and three-point loading, characterised by parameters a and b. From left to right the distance a of loading 
from the crack tip increases from 3 to 100. From up to bottom the asymmetry b of loading increases from to 0.99a. 



fl — 3 fl — 6 fl — 10 a = 100 




(St »t •f\ Kh 

Figure 5 : Shielding (light grey) , amplification (medium grey) and neutral (dark grey) regions created by a microcrack and 
a rigid line inclusion arranged as shown in Fig. [2^. The diagrams are for a bimaterial solid with rj = 0.67, and three-point 
loading, characterised by parameters a and b. From left to right the distance a of loading from the crack tip increases from 
3 to 100. From up to bottom the asymmetry b of loading increases from to 0.99a. 



for any symmetrical loading, independently of the distance a from the crack tip, and the diagrams are 
very similar regardless of the bimaterial parameter 77. For this reason, we present in Fig. [7] results only 
for the case of asymmetrical loading, 6 7^ 0, and weak interface between two identical materials, 77 = 0. 
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1*1 lf\ tfl Ifl 

Figure 6: Shielding (light grey), amplification (medium grey) and neutral (dark grey) regions created by a microcrack 
and a rigid line inclusion arranged as shown in Fig. [2^. The diagrams are for a bimaterial solid with 77 = —0.67, and 
three-point loading, characterised by parameters a and b. From left to right the distance a of loading from the crack tip 
increases from 3 to 100. From up to bottom the asymmetry b of loading increases from to 0.99a. 




<fl Vl Vl ¥>! 

Fi gure 7'. Shielding (light grey), amplification (medium grey) and neutral (dark grey) regions created by a microcrack 
and a rigid line inclusion arranged as shown in Fig. [2J3. The diagrams are for a weak interface between two identical 
materials, r] = 0, and three-point loading, characterised by parameters a and b. From left to right the distance a of loading 
from the crack tip increases from 3 to 100. From up to bottom the asymmetry b of loading increases from 0.33a to 0.99a.. 



It follows from Fig. [7] that the neutrality property of arrangement (b) is broken by the asymmetry of 
loading and thus amplification-shielding regions appear for b ^ 0. Moreover, the influence of asymmetry 
is not only quantitative but also qualitative, since the diagrams are completely different for different 
values of b. This emphasises the importance of the skew-symmetric weight functions (Piccolroaz et al., 
2009) which are used to evaluate the contribution of the skew-symmetric load. 
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4.2 Crack propagation and arrest 



In this section we discuss the influence of small defects on the crack propagation, in particular on the 
crack "acceleration" and arrest. Given a configuration of defects and position of the crack tip with 
respect to the defects, the formula ([5^ allows for the computation of the incremental crack advance 0. 
It is possible then to update the configuration with the new position of the crack tip with respect to the 
defects and recompute the incremental crack advance in the new configuration, following an iterative 
procedure. The crack "accelerates" when the increment (f) is increasing and "decelerates" in the opposite 
case. The total crack elongation is computed as 



where N is the number of iterations. 

In Fig. [5] we show the numerical results for the case of two defects (a microcrack and a rigid line 
inclusion) arranged as shown in Fig. [5^. The body is assumed homogeneous rj — and the loading is 
symmetric b ~ 0. The initial position of the microcrack relative to the main crack tip is characterised by 
the values c?i — I, ipi ~ tt/8. Five cases of angular orientation of the microcrack are considered, namely 
ai = 0, 7r/8, 7r/4, 37r/8, 7r/2, and denoted by labels from 1 to 5, see Fig.[SJi. The initial position and the 
inclination of the rigid line inclusion is determined by conditions (j67p . 

Fig. [H^ shows the crack elongation x as a function of the number of iterations. For all cases, we may 
observe crack acceleration followed by a rapid deceleration and arrest, see Fig. [5Jd. The crack arrests 
when a neutral configuration is reached at that particular position of the crack tip.. The total elongation 
at arrest is different for different cases. In the case denoted by 5, the crack advance is initially very slow, 
because the initial configuration is close to the border of the neutral region (dark grey region in Fig. 
[8ji), but the crack reaches the largest elongation at arrest. The incremental crack advance (j) is plotted 
against the crack elongation x in Fig. [S}:. 

In Fig.[5]we consider the same situation as in Fig.[51 however now the initial position of the microcrack 
relative to the main crack tip is characterised by the values di = I, ipi = 7tt/&, and four angular 
orientation are considered, ai = Stt/S, tt/2, bn/S, 37r/4, see Fig.[5Jl. We may observe a crack acceleration 
followed by a slow deceleration until a steady-state propagation is reached. In fact, as the crack is 
propagating and the angle ipi is increasing (up to the maximum value tt), the corresponding configuration 
will always remain in the amplification region, see Fig. [5Ji. 

5 Conclusions 

In the present paper we show that the new weight functions for a bimaterial plane constructed in Pic- 
colroaz et al. (2009) together with the dipole matrix approach provide an efficient tool in analysis of 
perturbation problems arising from the interaction of an interfacial crack with small defects, includ- 
ing elastic inclusions, microcracks, rigid line inclusions, and line defects with imperfect bonding. We 
also show that the skew-symmetric weight functions play a crucial role as they permit to evaluate the 
contribution of the skew-symmetric load applied along the crack faces. 

We illustrate the method with a number of applications, ranging from the analysis of shielding- 
amplification effects of the defects on the crack-tip fields to the perturbation modelling of the quasi-static 
interface crack propagation. The infiuence of asymmetry and position of the loading with respect to the 
crack tip is discussed and quantified. 

As a final remark, we note that the method can be extended to the case when the number of defects 
becomes increasingly high (clouds of defects). However, this requires additional accurate treatment of 
the interaction between defects and it is beyond the scope of the present paper. 



N 




(69) 
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Figure 8: Acceleration and arrest of the main crack along a weak interface between two identical materials, produced by 
a microcrack and a rigid line inclusion arranged as in Fig. [2ji. The initial position of the microcrack relative to the crack 
tip is characterised by the values di = 1, ipi = 7r/8. Five cases of angular orientation of the microcrack are considered, 
namely ai = 0, 7r/8, 7r/4, Stt/S, 7r/2, and denoted by labels from 1 to 5. The initial position and the inclination of the 
rigid line inclusion is determined by conditions I I67I I. The crack elongation, denoted by x, is assumed equal to in the 
initial configuration. Figure (a): crack elongation x vs. number of iterations. Figure (b): increment vs. number of 
iterations. Figure (c): increment vs. crack elongation x. Figure (d): initial configurations are indicated by spots in the 
shielding-amplification diagram. 




Figure 9: Initial acceleration followed by steady-state propagation of the main crack along a weak interface between two 
identical materials, produced by a microcrack and a rigid line inclusion arranged as in Fig. |2^. The initial position of the 
microcrack relative to the crack tip is characterised by the values di = 1, (pi = 7it/8. Four cases of angular orientation 
of the microcrack are considered, namely ai = Stt/S, 7r/2, Stt/S, 37r/4, and denoted by labels from 1 to 4. The initial 
position and the inclination of the rigid line inclusion is determined by conditions l|67|l . The crack elongation, denoted by 
X, is assumed equal to in the initial configuration. Figure (a): crack elongation x vs. number of iterations. Figure (b): 
increment <p vs. number of iterations. Figure (c): increment ip vs. crack elongation x. Figure (d): initial configurations are 
indicated by spots in the shielding— amplification diagram. 
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A APPENDIX 



A.l The dipole matrix for an elastic elliptic inclusion in antiplane shear 

Although the general problem of an elliptic inclusion in an infinite medium subject to remote uniform 
stress has been known since the mid of the last century (Donnell, 1941: Hardiman, 1954). the dipole 
matrix in the case of antiplane shear is practically impossible to find in the existing literature. For this 
reason, we have decided to present here not only the matrix but also in brief its derivation. 

We consider an elastic elliptic inclusion g, with shear modulus fj,Q, embedded in an infinite body with 
shear modulus fj,, undergoing antiplane shear deformation. The boundary of the inclusion is denoted by 
dg = {{x,y):xya^ + yyP = l}. 

The problem is formulated as follows: 

Alt = in \ 

Auq = in 5, 

du duo 
u = uo, = Mo^— on dg, 

on on 



u = Bbx - Bjy + 0{l) as r ^ ^/x"^ + y'^ ^ oo, (A.l) 

where uq and u are the displacements in the interior and exterior of the inclusion, respectively, n is an 
outward unit normal on dg and Bfj, Bj are known constants. 

The displacement in the inclusion is linear (Horgan, 1995; Ru and Scliiavone, 1996), so that 

Uq = Arx - Ajy in g. (A. 2) 

Our objective is to compute the next term in the asymptotic of the displacement field at infinity, eq. 
(jA.ip . which will be provide us with the required dipole matrix. To this purpose we will use the theory 
of analytic functions. Thus, we write uo{x,y) = KeU'^{z) = ^^{z) and u{x,y) = Re[/(z) = Ur{z), 
where U'^ = Az is a linear function in g and U is an analytic function in il\g satisfying the following 
conditions 

U ^ Bz + ^ + 0{\z\-^), B = BR + iBi, D = Dr + Wi, z ^ oo, 

UI-Ur^{), ^ioU°i-iiUi^{) on dg. (A.3) 

Note that, in order to find the dipole matrix, we do not need the complete solution but only the constant 
Di. Furthermore, we introduce the conformal mapping 

2 = c^C+^^, c> 0, 0<A<1, 

which maps the ellipse dg in the z-plane onto the unit circle in the C-plane. The segment d connecting 
the foci of the ellipse is mapped onto the circle |C| = "v/A, so that the region g\d\s transformed into the 
ring TA < Id < 1. 

We seek for the solution in the transformed coordinates in the form 

= Ac(C + A/C), C/ = /3C + 7i/C (A.4) 
Such a solution exists with constants j3 and 71 satisfying the following conditions 

_ ^i{l + x)-^iQ{l-\) Mi-A)-/io(i + A) 

^«^(l + A) + Mo(l-A)' ^'M(l-A) + /io(l + A)- 
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Going back to the original coordinates, from the asymptotic expansion C ~ z/c—cX/ z+0{\z\ '^), z — )• oo, 
we get 

/3 = cB, Di = 7ic - pXc. 



As a result, we obtain the required two terms asymptotic of the displacement at infinity in the form 

0, (A.5) 



u = Brx - Biy - ^{Br,, -Bi] ■ P + 0(r"'), as r 



27r ' a;-^ + y 

where, taking into account that c= (a + 6)/2 and A = (a — fe)/(a + 6), the dipole matrix P is given by 

/i - /io 



-Trab{a + h) 











A* - Mo 



(A.6) 



We obtain the solution we need in Sec. I3.1.f I by defining new functions u = u — Brx + Bjy and 
uq — uq — Brx + Biy. Then u and uo satisfy a BVP similar to (jA.ip . in which the condition at 
infinity is replaced by u — as r — > cx), and the ideal transmission conditions are replaced by imperfect 
transmission conditions as follows 



du 



diiQ 



u = uo, nj, Mo-^ 

on on 



(fi- ^0)^1 ■ {Br,- Bj} on dg. 



B Asymptotics at infinity for a movable rigid line inclusion 

One can find, as a particular case from Muskhelishvili (2008), a solution U for analytic fimction in the 
complex plane with given conditions along the symmetric segment of the length 2a situated along the 
real a;-axis: J/"*" = F{x), = F{x). The bounded solution for this problem reads 



U{z) 



^/r{z) r F{t)dt 



J-a y/mit-z)' 



(B.7) 



where z = x + iy, and R{z) = z"^ — . This solution vanishes at infinity if the following condition is 
satisfied 

r -4^=dt = 0. (B.8) 

It is clear from ([38)) that this condition is satisfied, since for a rigid line inclusion F(t) = Bt, where 
B = —'^xu^'^^Ys ' ^"^(0^3)^!; with $ being the rotation matrix and ei the unit vector along the x-axis. 
The asymptotics at infinity for this function takes the form 



Uiz) 



Ba^ 
~27 



+ Oiz~^), 



00. 



(B.9) 



Thus, the boundary layer solution for a rigid line inclusion in the coordinate system ^3 is given by a 
formula similar to ([27|l with A4i replaced by ([39]) . 
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